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Abstract 

Significant colour interference effects are expected in V{= 7, W, Z)+jet production 
at hadron colliders. These directly influence the hadronic antenna patterns and can 
provide a valuable diagnostic tool for probing the nature of the underlying parton 
subprocesses. Motivated by these ideas, we present quantitative predictions for the 
distributions of soft particles and jets in VF+jet production at the Tevatron pp collider. 



It was realised long ago that the overall structure of particle angular distributions in 
multijet events in hard scattering processes (the event portrait) is governed by the underlying 
colour dynamics at short distances. A natural idea has been proposed (for recent 

detailed studies see Ref. 0) to use this colour event portrait as a "partonometer" mapping 
the basic interaction short- distance process. The interest in this subject has recently been 
boosted because of two reasons. First, as was advocated in Ref. 0, the hadronic antenna 
pattern can be used as a diagnostic tool to dissect the colour structure of the exciting large- 
Ej- jet events observed by CDF |[T^ and DO |TT| at the Tevatron pp collider, as a way to 
distinguish between conventional QCD and possible new physics production mechanisms. 



Secondly, both CDF and DO have started very successfully [Il2|-p!5f| to measure the structure 
of multijet events and have demonstrated that the distinctive colour interference effects 
survive the hadronization stage and are clearly visible in the data. It has been known for a 
long time that an especially bright colour interference effect arises in the case of large- 
production of colour singlet objects, for instance, in V^+jet events (with V = 'j, or Z). 
The hadronic antenna patterns for to such processes are entirely analogous to that in the 
celebrated string |[T^ or drag [0 effect in e^e~ qqg events. 



Recently the first (very impressive) data on ly+jet production from DO [H, ^ have 
become available. The colour coherence effects are clearly seen and in our view these studies 
have a very promising future. They may play the same role for hadron colliders as the 
important series of results on inter-jet studies at e~^e~ colliders. 

The purpose of this paper is to present quantitative predictions for the colour interference 
phenomena in the distribution of soft particles and jets in V^+jet production at hadron 
colliders, in particular the Tevatron y/s = 1.8 TeV pp collider. 

There are two leading-order processes, qq — Vg and qg — > Vq. Each has its own 
distinctive antenna pattern, as we shall see. In principle, the antenna pattern could be used 
as a 'partonometer' to identify the dominant scattering process. 

There are two experimental methods which are likely to prove useful. With sufficiently 
high statistics, one could use an additional soft jet as the probe of the antenna pattern. 
Alternatively, with fewer events one could use the distribution of soft hadrons without re- 
quiring jet reconstruction. Both of these quantities are directly related to the inclusive soft 
gluon distribution in the NLO processes qq —>■ Vgg and qg — > Vqg. We will consider such 
distributions, first in the soft-gluon approximation where the results are particularly simple, 
and then using the exact QCD matrix elements and phase space constraints. 

The distribution of soft gluon radiation is controlled by the basic antenna pattern (see 
for example Ref. [0]) 

^ -p.-kprk- Elil - n, ■ n.) (1 - n, ■ n,) ' 

where the pf = Ei{l,ni) are the four-momenta of the partons participating in the hard 
scattering process. This corresponds to the emission of soft primary gluons with energies 
Ek- Aqcd ^ Ek <^ Ei,Ej. The particle fiow is then described by extra multiplicative 
cascading factors, see below. Since the parton scattering process acts as a colour antenna, 
the distribution of soft particle fiows between the jets is determined by the overall colour 
structure of the event [|l1-|^,p!^ 
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In order to make the study more quantitative we must define appropriate kinematic 
distributions and then compute the contributions of the various subprocesses. In order to 
fully understand the differences between these, we first consider the soft gluon distribution 
at the parton-parton scattering level with fixed kinematics. The generic process is 

a{pi) + h{p2) ^ c(p3) + d{pi) + g{k) , (2) 

where the gluon is soft relative to the two large- E't partons c and d. We assume massless 
quark and gluon partons and write My for the mass of the vector boson. In our quantitative 
studies, we are particularly interested in the angular distribution of the soft gluon around 
one of the large- ii^^ final-state particles, p^ say. Using the notation p^ = {E,px,py,pz), we 
write 

P3 = {Mt cosh ri,0, Et, MTsinh 7]) , 

= (A;j'Cosh(r7 + A?7), /cysin A0, fcycos A0, /crsinh(r7 + Ar])) , (3) 



where My = Et for c = q,g and Mt = y + My for c = V. The phase-space separation 
between the soft gluon and parton c is parametrized by At] and A0. Alternative variables, 
more suited to the experimental analysis, are the radial and polar angle variables in the 
"LEGO plot": 

A?7 = Ai?cos^, 

A(f) = Ai?sin/3. (4) 

which are defined in such a way that the LEGO-plot separation between soft jet k and hard 
jet/vector boson ^3 is R{k,p-i) = At]"^ + A02 = A_R, (0 < AR < 00), and the azimuthal 
orientation of k around p^ in the LEGO plot is parametrized by the angle P, {0 < f3 < 2tt). 
Our convention is such that (3 = corresponds to the direction of the incoming parton a, 
i.e. pi^ > 0. 

In terms of these variables, the soft gluon phase space is 

^ ^ krdkr ARdAR d(3 . (5) 



(27r)3 2Ek 167r3 

We will be particularly interested in the behaviour of the cross section as a function of /3 for 
fixed kT, AR and fixed Et, t]. 

The matrix elements for the lowest-order processes (for the case V = W) are:Q 



J2\M\\qipi)qip2) ^Wips)gip4)) 



gjg^ ( 1] t' + u' + 2sMl, 
A \ N^J tu 

-,2^2 1 +2 I ,2 I o„,/\#2 



j:\Mnq{pMp.)^w{pMP.)) = +s+2uM^ ^ 



^CKM factors are omitted for clarity. 
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where s = {pi +P2Y, t = (pi —psY and u = {pi —p^Y. In the soft-gluon approximation, the 
corresponding 2 — > 3 matrix elements are 

Y.\Mnqq - Wg{g)) = g^N^ ([U] + [24] - ^[12]j T.\M\\qq - Wg) 

Y^lM^qg - Wq{g)) = g^N^ [[12] + [24] - ^[14]^ T.\M\\qg - Wq) (7) 

Note that for these processes, the effect of the soft gluon emission is simply to multiply 
the lowest-order matrix elements squared by an overall factor consisting of three different 
antennae, defined according to (p, one of which is suppressed in the large N^. limit. This 
structure is universal for any electroweak boson + jet production, i.e. V = VF^,Z°,7,7*. 
Although we shall only study the V = W case in detail, similar results can be obtained for 
the other welLQ 

To illustrate the colour coherence properties of the above scattering processes, we first 
show results for a very simple kinematic configuration, 2 — >■ 2 scattering at 90° in the parton 
centre-of-mass frame, i.e. 77 = in Eq. @. For the case when the 'target' final-state particle 
(with momentum ^3) is a boson, the matrix elements squared are as given above and the 
various antennae which contribute to the radiation pattern are, in terms of the variables AR 
and P, 



[12] 
[14] 

[24] 



exp(A_Rcos/3) 



cosh(A/? cos (3) + cos(Ai? sin /3) 

exp(— Ai? cos P) 
cosh(A/? cos (3) + cos(Ai? sin (3) 



where an overall factor ofl/k^ has been omitted from the right-hand sides. When the target 
final-state particle is the quark or gluon, i.e. for the subprocesses qq — > gW and qg — > qW, 
the corresponding results are obtained by interchanging the momenta p^ and p^ in Eq. (|^. 

Fig. |l| shows the angular distribution of the soft gluon jet around each of the two 
large- £'t final-state particles, for AR = 1. The quantity plotted is 

^ ,2-2 I^I-^3P 



where the |A1nP are the 2 — n matrix elements squared of Eqs. and (^. Also shown for 
reference in Fig. |l] is the radiation pattern for the process qq VV, for which TZ = ACp = 
16/3, independent of j3. 

For the processes qq Wg and qq — >• gW, the distribution is symmetric about /? = 90° 
and much larger for the latter, as expected. Notice that the soft gluon distribution around 
the W in qq ^ Wg is smaller than that around the W in qq ^ WW; the away-side gluon 

^In particular, the results for Z°+jet production will be almost identical to those for M^+jet production, 
but in practice the event rates are an order of magnitude smaller. 
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jet in the former 'drags' the soft gluon away from the W hemisphere. In the formal hmit 
AR — >^ we find, for this simple kinematic configuration, 



n{qq ^ Wg) 




7 
3 



nQQ^gW) ~ (10) 

The divergent behaviour in the latter case reflects the matrix-element singularity when the 
final-state soft and hard gluons are coUinear. 

The processes qg — > Wq and qg qW have asymmetric f3 distributions in Fig. In 
both cases the distribution is larger in the backward {(3 > 90°) hemisphere, i.e. when the soft 
gluon moves in the same direction as the incoming gluon. In fact the distribution for the 
qg — >■ qW case is entirely analogous to that for e~^e~ — qqg, where the colour flow dynamics 
(string |]T6| or drag [|I^ effect) leads to a suppression of soft radiation between the outgoing 



quark and antiquark. 

The corresponding AR limits are (cf. Eq. (0)) 



n{qg ^ Wq) 




1 \ 22 



2m / 3 



nig^QW) ~ (11) 

We next address the question of whether the differences between the qq and qg subpro- 
cesses shown in Fig. ^ survive a more rigorous calculational treatment. We do this in two 
stages: first we compare the simple analytic soft-gluon approximation to the 2^3 ma- 
trix elements with the exact result, and then we include the effects of phase space, realistic 
experimental cuts, parton distributions etc. 

Fig. 1^ compares the exact and approximate calculations of the distributions of Fig. |T[ 
The exact 2 — > 3 matrix elements squared are calculated using the spinor techniques of 
Ref. |19|. Exact 2^3 kinematics are used, i.e. having fixed ps and k we define pi = —p^ — k 



and El = |p4p + m|. Results are displayed using three different values of the soft gluon 
transverse momentum kjn. Because the ratio of matrix elements squared is multiplied by k^, 
the analytic soft-gluon approximation (solid lines) is approached as kx — >■ 0. Evidently the 
soft-gluon approximation gives a very good representation of the exact result even at quite 
large values of kx/ET- 

The next step is to implement a full cross section calculation, including realistic cuts on 
the final state particles. We are interested in a final state consisting of a boson and two 
hadronic jets. The W and one jet are produced centrally with large Et, and the second jet 
is soft and required to lie in an annulus in the [t], 0) plane around either the W or the large- 
Et jet. The two jets are labelled according to the ordering of their transverse momenta, 
Et{J) > ExiJ). For illustration, we use the following set of cuts: 

Pt{W),Et{J) > 100 GeV , \r]{W)\, \r]{J)\ < 0.5 , 

Et{3) > 10 GeV , 0.7 < R{J or W,j) < 1.3 . (12) 
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Figure 1: Ratio TZ of the 2 — > 3 and 2 — »• 2 matrix elements squared times k^/g^ as a function 
of the soft gluon azimuthal angle about the large- jet or W, for the two subprocesses. 



In addition, we use the MRS(R2) (as(M|) = 0.120) parton distributions from Ref. [|T|]. One 
of the strong couplings is evaluated at scale /i = Et{J) and the other at /i = ExiJ). 

Fig. shows the (3 distributions calculated using the above cuts. The separate histograms 
correspond to the different initial states {qq and qg), and to the target particle (jet or W) at 
the centre of the annulus. The curves are normalized to the corresponding + 1 jet lowest- 
order cross sections for the corresponding initial state, with the cuts of ( |12D imposed. It is 
instructive to compare Figs. |^ and |^. The most noticeable difference is the suppression of the 
distributions at /3 = 0°, 180° in the full calculation. This has a simple kinematic explanation. 
For fixed values of the rapidities and transverse momenta of the W and large-i?T jet, the 
subprocess energy ^/I is maximal when the soft jet rapidity is largest, i.e. at /3 = 0°, 180°. 
These large \/§ configurations receive an additional suppression from the parton distribution 
functions. More importantly, the relative size of the various distributions is similar to that 
specified by the soft-gluon antennae and shown in Fig. ^ This is illustrated in Fig. ^, which 
shows (solid lines) the ratio of the 'jet' to W distributions in Fig. |^ for the qq and qg initial 




qq^WW : 



- qq^Wg - 
' ... I ... I ... I ... I ... I ... I ... I ... I ... ' 
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Figure 2: Ratio TZ of the 2 — > 3 and 2 — > 2 matrix elements squared times k'^/g'^ as a 
function of the soft gluon azimuthal angle about the large- jet or W, comparing the exact 
matrix element result (data points) for various values of the soft gluon with the fc^ — >■ 
approximations of Fig. |l] (solid lines). 



states. 

For the qg processes, the agreement between the full calculation and the soft approxi- 
mation is excellent]^ For the qq processes, the soft approximation appears to overestimate 
the result of the full calculation. However, this is due to the fact that the latter includes 
contributions from the 'four-quark' processes qq —>■ Wqq. Although these are formally sub- 
leading in the soft-jet limit [kx 0), they make a non-negligible contribution for this choice 
of parameters and cuts. Including only the qq Wgg processes gives instead the dot- 
ted histogram, and agreement with the soft approximation is restored. Fig. § shows the f3 
distribution when all subprocesses are included with their correct parton distribution weight- 

•^Note that in the full calculation AR is integrated over the range 0.7 < AR < 1.3 and the jet and 
W rapidities are integrated over the range \ri\ < 0.5, whereas in the soft approximation the antennae are 
calculated with AR — 1 and 77 = 0. 
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Figure 3: /? distributions for the + 2 jet cross sections at ^/s = 1.8 TeV, using exact 2^3 
matrix elements with cuts as described in the text, for different initial states and target final 
state particles. 



ing.0 For this choice of cuts, the qq and qg contributions are roughly equal, and so the net 
distribution interpolates between the separate distributions shown in Fig. ^ Note that in 
Fig. ^ the effect of including both qg and gq contributions is to symmetrize the distribution 
about P = 90°. Fig. ^ shows the corresponding prediction for the ratio of the jet and W 
distributions, i.e. the analogue of Fig. ^ but including all parton subprocesses and cuts. 

One problem with the analysis described above is that the event rate is very small. The 
requirement of a large Et central jet and a large pt central W strongly suppresses the cross 
section. Of course one solution would be to use a less restrictive set of cuts. Although this 
means that the asymptotic soft-gluon results are less applicable, the qualitative features of 
the colour dynamics are still apparent. As an example. Figs. ^ and || also show predictions 

"^The process gg — > Wqq is also included but makes a negligible contribution. 
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Figure 4: The ratios of jet/iy distributions from Fig. |]for the qq and qg initial states. The 
dashed hues correspond to the soft-gluon approximation, and the sohd histograms correspond 
to the exact calculation. In the dotted histogram, the qq —>■ Wqq contribution is absent. 



for a looser set of cuts: 

PT{W),ET{J)>30GeV , \7]{W)\, \7]{J)\ < 1.0, 

EtU) > 10 GeV , 0.7 < R{J or W,j) < 1.3 . (13) 

The effect of loosening the cuts is to slightly increase the relative contribution of the qq 
subprocess, which has the effect of enhancing the jet/PF ratio. The following table lists the 
values of this ratio for jets produced in the transverse plane (/? = vr/2) and in the event 
plane (/? = 0, tt). 





(3 = {0,n) 


(3 = n/2 


tight cuts 
loose cuts 


2.7 
3.2 


2.1 
2.4 
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Figure 5: /? distributions for the + 2 jet cross sections at ^/s = 1.8 TeV, with cuts as 
described in the text and including all subprocesses. 



A higher event rate could also be achieved by studying the distribution of soft jets in 
large Et 7+jet production, which should exhibit the same qualitative features as ly+jet 
production, although there may be experimental problems in obtaining a truly isolated large 
Erp photon sample. In fact, an additional theoretical advantage of 7+jet production is that 
there is a stronger dependence on Et of the subprocess decomposition, i.e. the relative 
amounts of qq and qg scattering. This arises because the average momentum fraction of 
the initial-state partons is (x) ~ 2ET/y/s, in contrast to (x) ~ {^J M^^r + E"^ + Et) / \/s for 
ly+jet production. Note that this argumentation only applies to central l^+jet production. 
In principle one can also study large Et events where both the V and the jet have large and 
comparable rapidities, such that one of the momentum fractions becomes large and the other 
small. The effect of this is to enrich the contribution from qg scattering, thus modifying the 
P distribution according to Fig. |[ 

An alternative approach, which requires much less luminosity, is to study the structure of 
particle (rather than jet) flows in the hard scattering events, see Refs. ||l|-|3|,p!7| . Here, unlike 
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Figure 6: The ratios of jet/iy distributions from Fig. 
the text. 



I for the two sets of cuts described in 



in the case of jet production, the geometry of the primary parton system is practically 
undisturbed, and there is of course no dependence on the jet-finding algorithm. Using soft 
particles, the colour coherence phenomena can be studied even when the overall statistics 
are limited. However, there is a price to pay. The success of the application of the analytical 
perturbative description to the particle flow distributions is based on the hypothesis of 
Local Parton Hadron Duality (LPHD) (see Ref. |]2D|), which assumes that the production 
of hadrons is governed by the QCD radiophysics of colour flows. The LPHD hypothesis is 
very well confirmed by all existing data, thus demonstrating that colour coherence effects 
successfully survive the hadronization stage, see Ref. pTl. Care must also be taken when 



using soft particles to study the colour coherence properties of the hard scattering to avoid a 
significant background contamination from soft particles produced in the 'underlying event'. 
In practice, this could necessitate a minimum px cut on the registered hadrons. 

Within the LPHD approach the radiation pattern for wide-angle soft hadron {h) produc- 
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tion can be written in the parton-parton centre-of-mass system as (see Refs. ^7^] 



where 



(14) 



(15) 



d\n Et 

and Ng^Ex) is the muhiphcity of hadrons h in an individual gluon jet with hardness Et- For 
the case when there is an experimental lower limit on the transverse energy of the registered 
hadron, 

Et > E:^^ > Aqcd , (16) 
then the last factor in Eq. ( |I4D becomes 



N'{E 



T 



h I 77imin\ 



NHE, 



(17) 



The cascading factor N^IEt) takes into account the fact that the registered hadron h is 
a part of the cascade initiated by a soft gluon jet g. It originates in the evolution equation 
for jet multiplicity to order ^/a^ (for details see Refs. 0, [T7[|): 



EdE. 



E. 



UN^^^N,{E,) 



27r 



Note that asymptotically 



N'(E) = -f N'(E) 



Nr 



with 



Ng{E) V 27r 



O 



(18) 

(19) 
(20) 



For instance, the exact result for charged particle production is (for details see Refs. I^O], ^ 

^ch 



iVf (E) = K^^ T{B) ( - 



-B+l 



Ib+i{.z) 



(21) 



with 



'16iV, 



r , F = In 



E 



Qo 



B 
b 



b I 3 ^ ^ 3iV2 



llNr - In 



f 



(22) 



Here I^, is the Modified Bessel function. The values of the parameters K^^ and Qq are 
obtained from fits to the e+e" data |21[] : 



K""^ ~ 1.3 , Qo - 250 MeV 



(23) 
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Finally, the distributions of charged particle flow accompanying the hard scattering processes 
qq Wg and qg Wq are given by 



Sirdnh 



Sndnfi 



1 r~~,\ dN, 



ch 
9 



dn 



n 



qq^Wg 



qg^Wq 



J141 + 1241-^112|J . 

[12] + [24]- ^[14] j (24) 



where 



(l-n-rii) (1-n-n,) ^ ' 

Note that these particle distributions take the form of an overall energy- dependent factor 
multiplying the same angular distributions as in the soft-gluon antenna patterns, Eq. (|^) 
and Fig. |l|. In studying the relative amounts of hadronic radiation near, say, the hard jet and 
the W, the dN^/dY factors cancel and the relative angular distribution can be obtained 
from the curves in Fig. |I|, weighted according to the relative contributions of the qq Wg 
and qg Wq subprocesses to the total ly+jet cross section. Fig. |^ shows the jet/VT 
/3-distribution ratio (i.e. the analogue of Fig. ||) for soft particle production. We assume 
here that the particles have much smaller transverse momenta than the hard jet and W, 
so that (i) there is essentially no phase-space suppression of the distributions, and (ii) the 
qq — > Wqq process makes a negligible contribution. In addition, the distribution is integrated 
over the annuliQ 0.7 < AR < 1.3 around the jet and the W, which are produced centrally, 
|?7| < 0.5. The curves correspond to different values of the hard jet Et- In fact for Ej- values 
in the range 20 GeV < Et ^ 100 GeV the qq and qg subprocesses give approximately equal 
contributions. The dashed lines in Fig. ^ are the limiting 'pure' qq and qg ratios. Notice 
that the particle jet/VT ratios are very similar in magnitude to the soft-jet ratios shown in 
Fig. ^. i.e. when expressed in terms of this ratio, the colour coherence effects appear to be 
universal. However, there is an important caveat to this conclusion. In this study we are only 
calculating the perturbative contribution to the particle and jet flow, i.e. the contribution 
associated with gluon emission from the single hard scatter which produces the W and the 
jet. In practice, there will also be a (non-perturbative) contribution from the 'underlying 
event'. The lower the particle/jet threshold, the more important this contribution is likely 
to be. We would further expect this background contribution to be approximately uniform 
in the [r], 0) plane, thereby increasing the (3 distributions by a constant amount. This would 
have the effect of decreasing the iet/W ratios shown in Figs. ^and|^. Since we are assuming 
that in practice the transverse momentum threshold for particle production will be much 
lower than for soft jet production, we would expect to see a larger effect in the former jet/lV 
ratio. 

In conclusion, the aim of this paper has been to exemplify the application of the radio- 
physics of colour flows to the partonometry of V{= 7, PF, Z)-|-jet events in hadron-hadron 
collisions. As was appreciated long ago [jl], 0], the colour interference effects accompanying 
such processes can be especially spectacular. They may well play for hadron colliders the 



^Recall that the curves in Fig. |i| are evaluated at AR = 1. 
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Figure 7: The ratios of jet/VT particle distributions, for different values of the jet = 
20, 60, 100 GeV. The dashed lines are the 'pure' qq and qg ratios. 



same role as the celebrated string [|1^] and drag [0 phenomena observed in three jet pro- 
duction in e~^e~ annihilation. Clear experimental observation of the effects discussed in this 
paper, and the quantitative agreement with the basic formulae presented here, would lend 
strong support to the idea that hadronic antenna patterns can provide a valuable diagnostic 
tool for elucidating the underlying dynamics in multi-jet high Et events. The first experi- 
mental results on PF+jet production from the DO collaboration |1^, ^ look very promising. 
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